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Abstract 

It's been recently demonstrated that quantum walks on graphs can 
solve certain computational problems faster than any classical algorithm. 
Therefore it is desirable to quantify those purely combinatorial proper- 
ties of graphs which quantum walks take advantage of and try and sep- 
arate them from those properties due to the encoding of the problem. 
In this paper we isolate the combinatorial property responsible (at least 
in part) for the computational speedups recently observed. We find that 
continuous-time quantum walks can exploit the covering space property 
of certain graphs. We formalise the notion of graph covering spaces. 
Then we demonstrate that a quantum walk on a graph Y which covers 
a smaller graph X can be equivalent to a quantum walk on the smaller 
graph X. This equivalence occurs only when the walk begins on cer- 
tain initial states, fibre- constant states, which respect the graph covering 
space structure. We illustrate these observations with walks on Cayley 
graphs; we show that walks on fibre-constant initial states for Cayley 
graphs are equivalent to walks on the induced Schreier graph. We also 
consider the problem of constructing efficient gate sequences simulating 
the time evolution of a continuous-time quantum walk. We argue that if 

Y > Ajv > Ajv-i > ■ ■ ■ — > Xi is a tower of graph covering 

spaces satisfying certain uniformity and growth conditions then there ex- 
ists an efficient quantum gate sequence simulating the walk. For the case 
of the walk on the m-torus graph T m on 2 n vertices we construct a gate 
sequence which uses 0(poly(n)) gates which is independent of the time t 
the walk is simulated for (and so the sequence can simulate the walk for 
exponential times). We argue that there exists a wide class of nontrivial 
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operators based on quantum walks on graphs which can be measured ef- 
ficiently using phase estimation. Interestingly, measuring these operators 
won't be unitarily equivalent to the quantum fourier transform. Finally, 
motivated by our results we introduce a new general class of computa- 
tional problems, HiddenCover, which includes a variant of the general hid- 
den subgroup problem as a subclass. We argue that quantum computers 
ought to be able to utilise covering space structures to efficiently solve 
problems from HiddenCover. 

1 Introduction 

There is a growing belief that quantum computers can solve certain computa- 
tional problems exponentially faster than any classical computer. Strong evi- 
dence for this belief comes in the form of Shor's algorithm Sho94 for factori- 
sation, and the graph traversal algorithm of Childs et. al. |CCD + 03) . 

Despite the spectacular success of the known quantum algorithms we believe 
that there is still only a fairly rudimentary understanding of the properties of 
quantum mechanics which are useful for computational speedups. We ascribe 
this poor understanding to at least two causes: (i) we are unsure what sorts of 
problems might be amenable to quantum-computational speed up; and (ii) even 
if we firmly believed that an efficient quantum algorithm existed for a problem, 
it is hard to come up with this putative algorithm because it is difficult (at 
least for us) to reason within the traditional quantum computing model — the 
quantum circuit model. 

All quantum algorithms are traditionally expressed in the quantum circuit 
model (see |NC00| and |Pre98j for a detailed description of the quantum circuit 
model). The quantum circuit model assigns a unit cost to certain elementary 
quantum gates, such as, for example, CNOT, H, and T (the 7r/8 phase gate). 
This is by no means the only way to express quantum algorithms. Three other 
computing paradigms are polynomially equivalent to the quantum circuit model: 
the quantum Turing machine model |BV97| . the one-way quantum computer 
RB01 , and the adiabatic evolution model (which has recen tly been sh own to 
be polynomially equivalent to the quantum circuit model in |AvDK + 04| '). 

Why might one want to consider computational paradigms other than the 
quantum circuit model? The reason is that there are conceptual peculiarities 
with the quantum circuit model which make it hard to work with when designing 
algorithms the traditional way. Firstly, there is a back-action that quantum 
superposition induces in gates like the CNOT where there can be a "backwards" 
flow of quantum information for certain initial states. Secondly, the existence 
of quantum entanglement seriously confuses the causal structure of quantum 
circuits because we can think of a Bell pair as sending a qubit backwards in time! 
However, in favour of the quantum circuit model is the fact that exponentially 
many degrees of freedom can be summarised succinctly. 

Some of the negative features of the quantum circuit model are amply ad- 
dressed in the quantum walk model of computation. Quantum walks, or the 
quantum dynamics of a free particle hopping on a graph, are an exciting new 
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paradigm for quantum computing. (For a review of quantum walks see |Kem03a| 
and references therein.) The attractiveness of quantum walks is that they pro- 
vide an extremely intuitive way to visualise exponentially many quantum de- 
grees of freedom. In contrast to the quantum circuit model, in a quantum walk 
there is a clear physical picture of where information is flowing and a definite 
notion of cause and effect. The price we pay for these intuitive features is that 
the exponentially (2™ for n qubits) many degrees of freedom of n qubits trans- 
late to exponentially vertices in the graph. Additionally, there is no clear way 
to translate a quantum walk efficiently into the quantum circuit model. Despite 
these difficulties we believe that quantum walks provide an attractive method- 
ology in the quantum algorithm designer's toolkit because they appeal to the 
geometrical intuitions. 

The growing number of quantum-walk algorithms might be seen as evidence 
for the conceptual utility of the quantum walk model. We mention three re- 
cent algorithms: (i) the quantum walk search algorithm SKWQ H], ICG03 , 
and |AKR04j : (ii) the graph-traversal algorithm of Childs et. al. |CCD + 03 ; 
and (iii) the element distinctness algorithm of Ambainis Amb03 , and relatives 

|HEn^lM^n3ll^04] . 

What unites the quantum walk algorithms? (And, more ambitiously, all 
quantum algorithms?) Let's concentrate on the results of |CFG02I lCCD + 03l 
MR02 , Kcm03b , which are based on the continuous-time quantum walk. (Note 
that the quantum walks on the hypercubc MR02, Kem03b and the graphs in 
CFG02 don't provide the speedups for the solution to any algorithmic prob- 
lem.) The results in these papers appear to be related phenomena (they all 
take advantage of "column spaces"), however, this relationship has not yet, to 
the best of our knowledge, been quantified fully. In this paper we identify and 
generalise the combinatorial property of these graphs which leads to small hit- 
ting times. We believe this is a key ingredient underlying quantum speedups of 
hitting times. The combinatorial property we isolate is that all of the graphs 
walked on in these papers are covering spaces for much smaller graphs. 

Quantum computers and covering spaces 

What is a covering space? Suppose that X and Y are arcwise-connected and 
locally arcwise-connected topological spaces, respectively. Then (Y, tt) is said to 
be a covering space of X if tt : Y — > X is a surjective continuous map with every 
x G X having an open neighborhood U such that every connected component 
of 7r _1 (i7) is mapped homeomorphically onto U by tt. The preimage of a point 
in X is called a fibre of tt. An example covering space is shown in figure ^ 

In this paper we will show that a continuous-time quantum walk on a graph 
Y which is a covering space for another graph X (in the natural topology) 
is isomorphic to a quantum walk on X, under certain specific circumstances. 
This equivalence occurs when the walk starts on fibre- constant quantum states. 
Additionally, we will argue, and show in some cases, that if there is a tower of 
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Figure 1: An example of a covering space. In this case the circle X = S 1 is 
covered by another circle Y which has twice the circumference of X under the 
projection n. Note that the inverse projection of any point x e X is a finite 
set consisting of two points a and b in the covering space Y. One can think of 
this covering space as a subset of the Riemann surface for f(z) = ^fz, which is 
a covering space for C. 



covering spaces 

Y > A N > Ajv-i > ■ ■ ■ > Ai [1) 

then there is an efficient gate sequence (i.e. using 0(poly(n)) elementary gates) 
for the quantisation of Y. The idea we exploit is to recursively and hierarchi- 
cally construct the gate sequence from the "elementary gate" U(Xi) and the 
specification of how U(Xi) lives in U(X2 )■ 

We believe that these results are not specific to continuous-time quantum 
walks, but rather indicative of a general principle. 

We formulate this principle in the following way: consider some collection 
of mathematical objects like the class of simple graphs, or something with more 
structure, like the class of finite fields. Let U(X) be a quantisation scheme for 
this class of objects, by which we mean a way to associate a unitary matrix 
U(X) acting on a finite-dimensional Hilbert space Ti-(X) with any object X. 
Suppose, further, that Y is a covering space 7r : Y — > X for another object 
X. If the quantisation scheme U(X) is "sufficiently well-behaved" then U(Y) 
should be related to U(X) according to a map, the pull-back ir* : U(X) — > 
U(Y). If, further, Y is determined from A" in a sufficiently simple way, then it 
should be possible to construct U(Y) from U(X) according to this specification. 
Mimicking the recursive and hierarchical construction for gate-sequences for 
continuous-time quantum walks we expect that if 

v TJV Y TN-l v 7TJV-2 TTl „ 

Y ► Ajv ► Aat_x ► • • • > Ai (Z) 

is a tower of covering spaces obeying some uniformity and growth conditions 
then there should be an efficient gate sequence for the quantisation U(X). 
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We retroactively argue that there is an example of where this philosophy has 
already proven successful. In this example we take as our space of mathematical 
objects the category GrpFin of finite groups. We take the quantisation scheme 
to be the contravariant functor which associates the quantum fourier transform 
with every finite group. (For an introduction to category theory see |ML98| .') For 
certain towers of subgroups, i.e. those which are polynomially uniform MRR03 , 
there is an efficient hierarchical scheme to construct the quantisation U (Y) . Here 
we are taking the "is a subgroup" relation to be the covering space relation in 
this category. To accord this exactly with the definition of covering space given 
earlier requires the introduction of an appropriate topology on a finite group. 
This is something we will avoid, preferring only this qualitative argument - 
clearly the analogy isn't perfect. 

What sort of problems could these supposed efficient quantisations solve? 
We will argue that they could be used to solve a problem we call the Hidden 
Covering Space problem which includes the hidden subgroup problem as a sub- 
class. Roughly speaking, if there is a function / on Y which is periodic on some 
object X covered by Y, then a quantum computer should be able to identify 
the structure X efficiently. 

The general properties of covering spaces are extremely tantalising and sug- 
gestive. We hope to convince the reader that the results we have found concern- 
ing graph covering spaces and continuous-time quantum walks are indicative of 
a much larger framework. We believe that many quantum algorithms exhibit- 
ing an exponential separation between classical and quantum complexity could 
exist. We argue that the properties of quantisations of covering spaces may 
provide many opportunities to design such new algorithms. This is not least 
because covering spaces are pervasive in mathematics, from number fields and 
algebraic surfaces to geometric group theory, and, of course, Galois theory. 

The outline of this paper is as follows. We begin in by describing the 
quantisation scheme we study in the remainder of this paper, the continuous- 
time quantum walk. In ^Jwe review the theory of graph covering spaces and the 
heat kernel for graphs. We then apply these results in jJl|to demonstrate that 
quantum walks which begin on fibre-constant states are isomorphic to quantum 
walks on smaller graphs. We illustrate our results in for the hypercube and 
Cayley graphs. In S^we utilise the covering-space properties of the m-torus 
graph on 2™ vertices to construct an efficient gate sequence for their quanti- 
sation. Motivated by our results we introduce, in a new class of problem, 
HiddenCover, which quantum computers may be able to solve efficiently. We 
also solve the hidden cover problem for the m-torus graph which, incidentally, 
provides a (philosophically) different way to solve the abelian hidden subgroup 
problem. 
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2 Quantisations of Graphs: the Continuous-Time 
Quantum Walk 



In this section we introduce the continuous-time quantum walk, which is a quan- 
tisation scheme for the class of simple graphs. For further details about the 
graph-theoretic notation and terminology we use in this section and the rest of 
this paper see |Big93| IHDSMl IChu97| . 

Let us begin by defining the main objects of our study. By a weighted 
graph Y we mean a vertex set V = V(Y) with an associated weight function 
w : V x V — > R + (M + denotes the real numbers i£l such that x > 0) which 
satisfies 

w(u, v) = w(v, u). (3) 

If w(u, v) > then we refer to {it, v} as an edge of Y, and we say that u and v 
are adjacent. By a simple graph we mean the special situation where w(u, v) is 
either or 1 and w(u, u) = for all u € V. 
We define the degree d v for a vertex v to be 

d v = ^2 w ( u , v )- (4) 

A graph is regular if all the degrees are the same. 

The Laplacian of a weighted graph Y on n vertices and weight function w 
is the n x n matrix A given by 

!d v — w(v, v) if u = v, 
—w(u, v) if u and v are adjacent, (5) 

otherwise, 

where u and v are two arbitrary vertices in V. 

We now consider introduce the (formal) complex vector space Ti.(Y) spanned 
by the vectors \u), u £ V, which we take to be orthonormal under the natural 
inner product: (u\v) = 5 UV . Of course this vector space is isomorphic to C™. 
We think of the vector space TLiY) as the space of all complex- valued functions 
/ on the finite set V(Y), where for each vector |/) = X)uev(Y) fu\ u ) G 
we define the function / : V(Y) — » C by f(u) — (u\f). 

We define the adjacency matrix of a weighted graph Y to be the operator 

A0r)= E H^v)\u)(v\. (6) 

u,v£V(Y) 

The Laplacian, acts in a natural way on "H{Y) as 

A= E {d u -w(u,u))\u)(u\- E w ( u > v )\ u )( v \> ( 7 ) 

uGV(Y) u,vEV(_Y) 
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where we use the notation u ~ v to mean that u is adjacent to v and u ^ v. 
For a specific vector |/) = J2 u ev(Y) fu\ u ) we have 

A|/)= Yl (m-f(v))w(u,v)\u). (8) 
u,«£V(y) 

It is worth noting that the Laplacian can be written A = D(Y) — A(Y), where 

D(Y) -I2u,veV(Y) W ( U ' V )\ U )( U \- 

We define the heat equation and Schrddinger equation for Y to be the dif- 
ferential equations 

J^(r))=-A|#r)>, (9) 

and 

i^m))=&\m), (io) 

respectively. Note that the Schrodinger equation is equivalent to the heat equa- 
tion with the replacement r = it. 

The Laplacian A for a graph on n vertices is a symmetric matrix and so we 
can write its spectral decomposition 

n-l 

A = ^A i |A J -)(A j |. (11) 

3=0 

We often refer to an eigenstate of A as a harmonic eigenfunction. 

The spectral decomposition of the Laplacian allows us to define the heat 
kernel and propagator 

n-l 

H(r) = Y / e-^\X j )(X j \, (12) 

3=0 



and 



ii-i 



C/W = ^e-^*|A J -)(A J -|, (13) 
i=o 

respectively. Note that U(t) — H{ir). Using the the heat kernel and propagator 
we can solve the heat and Schrodinger equations (j^J) and (|l(Jf> with initial state 
|V(0)) by defining |^(r)) = ff(r)|^(0)) and |^(t)) = C7(t)|V(0)>: 

(14) 

= -A#(r)|V(0)). 

The result for the propagator follows by substituting r = it. 

Definition 2.1. A continuous-time quantum walk on a weighted graph Y is 
the propagator U[Y](t) associated with the Laplacian A[Y] for Y. 
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3 Graph Covering Spaces and the Heat Kernel 



In this section we review the theory of graph covering spaces and the heat kernel 
for weighted graphs. For an introduction to graph covering spaces see Big93 
and |GM80j. For further results i n this area see |UY99| and [5T96ll5ToTllTer99l 
ISTT)01ITct02] . We follow |UY99| closely for most of this section. The principle 
result of this section is that the spectrum of a graph X covered by another 
graph Y is contained in the spectrum of Y. The eigenvectors of X also induce 
eigenvectors in Y . 

We begin with some definitions. 

Definition 3.1. Suppose we have two graphs X and Y with weight functions 
u>x(x, y) and wy{u, v), respectively. We say Y is a covering space 1 for X (alter- 
natively, X is covered by Y) if there is a set map ir : V(Y) — > V(X) satisfying 
the following two properties: 

1. There is a \i € R + (the positive real numbers including 0), called the index 
of 7r, such that for u, v E V(X) the equality holds: 

^ w Y (x,y) = fiy/\Tr- 1 (u)\\TV- 1 (v)\wx{u,v). (15) 

a;G7r _1 (u) 



2. For x, y G V(Y) with tt(x) = n(y) and v € V(X), we have 

^2 w Y (z,x)= ^2 w Y(z',y)- (16) 



We now translate Definition 13.11 into a more useful statement about the 
adjacency matrices of Y and X. In doing so we make contact with the theory 
of equitable partitions |Big93| . 

We begin by considering an arbitrary projection operator tt : V(Y) — » V(X) 
and define the corresponding pull-back operator P : H(Y) — > H(X), 

uev(x) Vk («)l 

(It is easy to establish the identity (P'P) 2 = P^P via multiplication, which 
shows that P^P is a projector. Similarly, one can show that PP' = Iu(x)- 
This implies that P* is an isometric imbedding pt : W(X) -> H(y).) What 
are the conditions that P must satisfy in order that 7r be a legal covering map? 
We answer this question in two steps, by translating parts (i) and (ii) of Defi- 
nition into conditions on P. 

1 Note that our definition of (unramified) graph covering spaces accords with the standard 
definition given in the introduction when we take for the topology of the graphs X and Y the 
weak topology. Sec Hat02 and ST01 STOO for further discussion and details. 
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Using the operator P we define the quotient matrix A 7r {X) = PA(Y)P\ 
which we want to identify with the adjacency matrix of X. Expanding this 
expression and utilising Definition 13.11 we obtain 

mx)= E E E /, Zt v \ ( ,, mm\m){v\ 

uev(x) vev(x) x , ye v{Y) VF WlF H 

-EE , WY } l ' m) Nfol (is) 

len 1 (u)m£ir 1 (v) 

= E E l*u>x(u,v)\u){v\= fiA(X). 

uGV(X) veV(X) 

This expression shows us that the adjacency matrix wy(u,v) for Y satisfies 
the condition (i) of Definition 13.11 to be a covering space for the graph X with 
adjacency matrix wx(u,v). Note that part (i) of Definition 13. II is really saying 
that Y is a covering space for any graph which has as its adjacency matrix a 
scalar multiple of PA(Y)P <! . 

In order to establish the matrix version of part (ii) of Definition 13 . II we need 
to make use of the following result of Godsil and McKay GM80 , which is the 
fundamental connection between P and graph covering spaces. We include a 
proof for completeness. 

Lemma 3.2 (Godsil and McKay [GM80 _). The graph Y is a covering space 
for X if and only if PA(Y)=A{X)P. 

Proof. For a vertex v £ V(Y), 1 < j < |V(X)| define h V j to be the sum of the 
weights of the edges connecting vertices in 7r -1 (itj) with v, where Uj € V(X). 
That is, h vj = E^-^ttfi («,«)■ For 1 <j < \V(X)\, v € V(Y), we have 

(uj\PA(Y)\v) = 1 M*,v)= , hv l , (19) 



and 

( Uj \A(X)P\v) = ^d. 



where v S n (ttfe). 

By comparing (JUJ and we note that if PA(Y) = A(X)P, then h VJ 

equals J j^Zi ^"^ jtt>2 (v,j , Uk) for all v E tt -1 ^-) and so ir satisfies part (ii) of 
definition 13.11 

For the converse, suppose it satisfies the conditions of Definition 13. II We see 
that 

PA(Y) = PA(Y)P^P = A(X)P. (21) 

□ 
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It is well known that the spectrum of a graph Y which is a covering space 
for X is determined, in part, by the spectrum of X. Let's make this a bit more 
precise. Suppose we have a harmonic eigenfunction of \f) of X with eigenvalue 
A. We can lift or pull back this eigenfunction to a harmonic eigenfunction \g) 
of Y, by defining for each vertex I £ Y, (l\g) = (7r(Z)|/)/\/7r _1 (Z). Equivalently, 

\g) = P i \f)- 

We use the notation \tt 1 (u)) = l/-\/|7r _1 (tt)| E/e-n— 1 (u) 10 anc ^ WT &e 

\g)= E /(0k _1 (0>. (22) 

leV(X) 

It follows from Definition 13. II that \g) is a harmonic eigenfunction for Y. 

Lemma 3.3. Let Y be a covering space for X with projection map it. Then for 
any eigenvector \ f) and scalar A, A(X)\f) = X\f) if and only if A(Y)P*\f) = 
AP f |/). 

Proof. If A{X)\f) = A|/), then P^A(X)\f) = APt|/) and so A{Y)P^\f) = 
AP+I/), by lemma EH If A(Y)P*\f) = AP+|/), then PA(F)Pt|/) = APP^]/) 
and so A(JC)|/) = A|/). □ 

We define the Laplacian A(X) for the covered graph X via A(X) = PD(Y)P^- 
PA(Y)Pt. 

The main result for this section is the following lemma relating the heat 
kernels and propagators of a graph X and an arbitrary covering space Y. 

Lemma 3.4. Suppose Y is a covering space for X. Let H[Y](t) and H[X](r) 
and U[Y](t) and U[X](t) denote the heat kernels and propagators for Y and X , 
respectively. Then we have 

H[X]{t) = PH[Y](t)P\ and U[X](t) = PU[Y](t)Pl (23) 

Proof. We establish the lemma for the heat kernel. The result for the propagator 
follows from substituting r = it. 

Note first that A(X) r = P A(Y) r P\ which follows from a simple induction 
using lemma EOl 

Expanding the heat kernel in an absolutely convergent power series in r 
gives: 



Z J n | 



(24) 

f^ptl^pLpi = p H[YK r)pi. 

3=0 J ' 

□ 
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Because we are living in finite Hilbert spaces, where everything is well- 
behaved, we don't have to worry about the convergence of series. For this 
reason, knowing the heat kernel H (r) is equivalent to knowing the propagator 
U (t) . We take advantage of this fact by only stating all our results about the 
propagator in terms of the heat kernel in imaginary time. 



4 Induced Quantum Walks — the Quotient Walk 

In this section we show that if a graph Y is a covering space for another graph 
X with projection tt : V(Y) — > V(X) then a quantum walk on Y which begins 
on a state constant on fibres of tt is isomorphic to a quantum walk on X . This 
result provides an insight into the h itting-tim e speedups observed by Kempe 
Kcm03b and Childs et. al. ICFG02I ICCD^ . 



Definition 4.1. Let Y be a covering space for X with projection tt : V(Y) — > 
V(X). We say that a quantum state \ip) S H(Y) is fibre- constant for tt if 
can be expressed as 

IV) = Y c «k _1 W). ( 25 ) 

uEV(X) 

where we are using the notation 

\n-\u))=P^\u)= 1 Y, !*> ( 2fi ) 

of the previous section. We write \ip) = P'\<f>), where \<f>) — Y^ueV(x) c u\ u )- 

What is the time evolution of a fibre-constant state \ip) for covering map tt? 
The following simple lemma shows us that it is isomorphic to a walk on X . 

Lemma 4.2. Suppose Y is a covering space for X with covering map tt. Let 
\tp) = P*\4>) be a fibre- constant state. Then the time evolution of\ip(t)) on Y is 
isomorphic to the time evolution of \(f>) on X . 

Proof. Let H[Y](t) be the heat kernel for Y. The time evolution of \ip(r)) on 
Y is given by 

\iP(r))=H[Y}(T)pi\<f>} 

= P^H[X] (r) | <f>) by lemma EH1 

This equation makes it clear that the time-evolution of \ip[r)) is determined 
from that of \<b) on X. □ 



5 Examples: the Hypercube Q n and Cayley and 
Schreier Graphs 

We now illustrate this result for the hypercube and the Cayley and Schreier 
graphs of a finite group. 
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Definition 5.1. Let G be a group, and let S C G be a set of group elements 
such that: (1) The identity element e ^ S; (2) If x G S, then a; -1 € S. The 
Cayley graph X(G, S) associated with G and S is then defined as the simple 
graph having one vertex associated with each group element and directed edges 
(<7, h) whenever gh^ 1 £ S. (It is easy to check that with our definition of the 
generating set this graph is well-defined.) 

Example 5.2. We think of the hypercube graph Q n of dimension n as a 
Cayley graph of the the abelian group (Z/2Z)™, the n-fold direct product of 
Z/2Z. Specifically, we write Q n = X((Z/2Z) n , {ei, e 2 , . . . , e„}), where ej = 
(0, . . . , 1, . . . , 0) is the unit vector with a one in the jth entry. Note that 
^(Qn)! = 2". It is easily verified that Q n is a covering for the weighted path 
of size (V'(P^)I =71+1 with adjacency matrix 



y/E 
v/n y/2(n-l) 

y/2(n-l) 











(28) 



The projection map ir : V(Q n ) — > V(P®) maps the collection of vertices of 
V(Q n ) of hamming weight j to the single vertex vj £ V(P$) at position j. 

Remark 5.3. Because Q n is vertex-transitive (i.e. the automorphism group 
Aut(Q„) acts transitively on V(Q n )) a quantum walk on Q n beginning on the 
state concentrated on any single vertex is equivalent to a walk on the weighted 
path pQ beginning at the first vertex. (Note that all Cayley graphs Y are 
vertex-transitive, where the automorphism group Aut(y) acts transitively on 
the vertex set.) The adjacency operator f2"%| of the weighted path is exactly 
the x angular momentum operator J x in the spin-^ irrep. of 517(2). One can 
envisage the quantum state propagating from the first vertex to the last vertex 
in P$ as a localised state at z = +| on north pole of the Bloch sphere prop- 
agating across the sphere (and hence becoming completely delocalised around 
the equator in the process) to the south pole z = — 

Our second example concerns the Cayley graph for an arbitrary finite abelian 
or nonabelian group G. To introduce it we need a definition. 

Definition 5.4. Given a subgroup H of a finite group G and a generating 
subset S (recall s S S implies s _1 € S and e £" S) of G, which we call the edge 
set, the Schreier graph X = X(G/H, S) has as vertices the cosets gH, g e G. 
Two vertices gH and s~ 1 gH are joined by an edge for all seS. 

Remark 5.5. Note that unless H is a normal subgroup H < G then the Schreier 
graph X(G/H, S) will contain loops. In the case that H is normal then X(G/H, S) 
is exactly the Cayley graph of the group G/H with generating set S. Note also 
that our definition of graph covering spaces includes graphs with loops (i.e. when 
w(u, u) > 0) so that the following results are well-defined. 
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Lemma 5.6. The adjacency matrix A{X) for a Cayley graph X(G,S) can be 
written in the following way 

A(x) = j2J2\ x ^ s ' lx \- ( 29 ) 

Similarly, the adjacency matrix A(X(G/H, S)) for the Schreier graph X(G/H, S) 
is given by 

A(X(G/H,S))=J2 E m^gEY (30) 

seS gHeH 

We aim to show that if Y is a Cayley graph Y = X(G, S) for a finite group 
G with generating set S then it is a covering space for the Schreier graphs 
X(G/H, S) for all H < G. In order to show this we define the projection map 
7r : Y — > X(G/H,S) by ir(g) = gH. The corresponding pull-back operator is 
then given by 

p = E -4m\9H)(g'\. (31) 



gHeG/H 

g'egH 

In order to establish our claim we need to show that PA{Y) and A(X(G/H, S))P 
are equal and apply lemma 

A{X{G/H,S))P=-^=y j V V \gH){s- l gH\kH){k'\ 



\E\ seS gHeG/H kHeG/H 
k'ekH 



t=E E *.-',wiff*><*'i 

I I ses kHeG/H C32) 
1 ^ ^ IfcHXa- 1 * 7 ! 



VWl s es kHec/H 

k'ekH 

= PA(Y). 

The transition from the second-last line to the last follows by a change of variable 
kH i— ► s&ff and fc' i— * sk', and the definition of P and A(Y). This discussion, 
together with lemma constitutes the following corollary. 

Corollary 5.7. A quantum walk on a Cayley graph X(G, S) which begins on the 
fibre- constant state \4>h) = y/\H\/\G\ ^2gHeG / h c gH\' K ~ 1 (gH)) is isomorphic to 
a quantum walk on induced the Schreier graph X(G/H, S) for the subgroup H 
which begins on the induced state \<j>') = y/\H\/\G\ ^2 g HeG/H c gn\gH) . 

6 Efficient Gate Sequences for Continuous-Time 
Quantum Walks 

In this section we explain how to construct efficient gate sequences simulating 
continuous-time quantum walks on certain classes of graphs Y. The principle 
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feature intuitively exploited in these gate sequences is that the graph Y covers 
a tower Y -~^-> Xj^ N 1 > X^_i N 2 > • • • — U X\ of graphs Xj. An interesting 
feature of these gate sequences is that their length doesn't depend on the time 
t the walk is simulated for. 

Consider the cycle on 2 n vertices, which is the Cayley graph = 
X(Z/2"Z, {±1}) of the cyclic group. We begin by presenting a quantum circuit 
which simulates the quantum walk on the cycle C2™ which uses only 0(poly(n)) 
gates. A special feature of this gate sequence is that it can simulate the walk 
for any given time t, including exponential times 0{t) = 2 poly (") . 

Recall |Lub95| that the eigenvalues and eigenstates of the adjacency matrix 
A(C m ) are given in terms of sums of the characters /i J = e l ~ 3 , j = 0, . . . , m — 1, 
of Z/mZ. Specifically, 

m — 1 

A{C m ) = X] W + ^ j )\W(j))(W(j)\, (33) 
j=o 

where 

1 171 — 1 

\W(j)) = ^J2^ k \k). (34) 

Note that the vectors | W(j)) are precisely those given by applying the quantum 
fourier transformation to \j): 

m— 1 

\W(j)) = QFT|i) = -f= M W |*>» (35) 



Our objective is to simulate, with a quantum circuit, the evolution 

m — 1 

U[C 2 n](t) = £ e- j2cos ^ J ) t |W r (j))(^(j)|, (36) 
i=o 

where m = 2™. We note that this unitary can be written as 

U[C 2n ]{t) = QFT<f>(t)QFT t , (37) 

where 

m— 1 

*(*) = E ^ 2coB ^>|i)(i|. (38) 

The phase operation $(£) can be performed efficiently using, for instance, the 
phase kickback trick CEMM98 . The phase kickback trick shows that if the 
computation 

\j) ~ mm), (39) 

where \j) is a computational basis state of n qubits and f(J) is an n bit ap- 
proximation to (2 cos (— i) i) mod 27r (calculable efficiently classically), is im- 
plementable efficiently, then the phase changing operation 

|j)~ e - iac »(^)*|.7) J (40) 
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is implementable using 0(poly(ra)) quantum gates JCEMM98 . (It is easy to 
obtain arbitrary accuracy by appending more qubits to the register for |/(a;)).) 

Hence, because there is an efficient (i.e., using 0(poly(n)) elementary gates) 
quantum circuit for the quantum fourier transform |NC00| . there is an efficient 
gate sequence approximating the propagator U[C2^](t). Moreover, the number 
of gates required to simulate the propagator does not depend on the time t. 

We now note some features of our gate sequence for U[C m ](t). Firstly, we 
point out that the circuits for U[C m ](t) can be easily generalised to simulate 
quantum walks on a many other families of graphs. This is because any circu- 
lant matrix is diagonalised by the discrete fourier transformation, and so has 
the same eigenvectors as C m . This observation allows us to generalise our sim- 
ulation sequence to any circulant matrix whose eigenvalue set {Xj} is efficiently 
calculable using classical or quantum means. In particular, this means that all 
Cay ley graphs of Z/mZ and the Paley graphs 2 are efficiently simulatable using 
relatives of the quantum circuit we have constructed. Additionally, any carte- 
sian product of graphs which are efficiently simulatable via the gate sequence 
above will be efficiently simulatable. This is because the adjacency matrix for 
the cartesian product of two graphs X and Y is A(XxY) = A(X)®I+I® A(Y). 
In particular, a gate sequence for the m-torus T%i on 2™ vertices follows from 

m — 1 

e -iA<r*)t = Q e -iA{C a n)t t (41) 
3=0 

The second feature we would like to point out is that the covering space 
properties of C^n, in particular that C%n is a covering space for C 2 n/2, allow us 
to "rediscover", using geometric constructions, an efficient quantum circuit for 
the quantum fourier transform. Roughly speaking, an eigenvector of the base 
graph C 2 n/2 can be lifted to 2™/ 2 eigenvectors of C2™ via a simple unitary recipe. 
This is a general feature of the laplacian on a covering space, and extends even 
to the continuous case |Ros97| . This idea is explored further in |Osb04j . 

We believe that the existence of covering spaces structures in graphs may be 
utilised more generally to supply efficient quantum circuits for quantum walks 
on graphs. We sketch one result which is indicative of this idea. 

We restrict our attention to regular graphs. Let Y be a regular graph on 
2 n vertices which is a covering space for a tower of graphs Y — X^ - 1 > 
X N -i ■■■ ^ X 1 with the property that \V{X ]+1 )\ w y/\V(Xj)}. In this 

case N must be 0(poly(log(|V(y|))). 

We assume that m = \ttJj (u)\ is the same for all u £ V(Xjy). Consider the 
transition matrix A{u,v) for an edge (u,v) € E(Xjy), which is the 2m x 2m 
matrix whose (x,y) entry, where x,y S 7r^ 1 (w) U 7r_^- 1 (?;), is given by A{Y) x , y . 

2 Given a finite field ¥ q with q elements, the Paley graph P(F q ) is the graph with vertex 
set V = ¥ q where two vertices are joined when their difference is a square in the field. This is 
an undirected graph when q is congruent l(mod4) IChu97l . Note that Paley graphs have edge 
density approximately i , and so are not efficiently simulatable using the recipe of [XTS03 . 
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Suppose we can label the vertices in tt n (it) so that 

A(u,v)=\? Qj®I m , or A(u,v)=(° l ^\<E)^=J m . (42) 



where I m is the m x m identity matrix and J m is the m x to all l's matrix. 
We call such a transition matrix trivial, and the we say that the edge (it, i;) has 
trivial pull-back. 

Suppose the edges of A{Xm) have trivial pull-back for all but 0(poly(log(| Y(Y)|))) 
edges. In the case we can write the adjacency matrix for A(Y) as 

A(Y) = A(X N ) I m + V, or A(Y) = A(X N ) <E> -^J m + V, (43) 

where T> is a symmetric |V(Y)| x |V(Y)| matrix whose (x,y) entry, x,y £ 
■njf (u) U 7r^ 1 (w), is nonzero only when there is an edge (u,v) g E(Xn) whose 
transition matrix is nontrivial. When the edges of A(Xn) have trivial pull-back 
for all but 0(poly(log(|V(Y)|))) edges then V is s parse, a nd has 0(poly(log(|V(Y)|))) 
nonzero entries in each row (in the language of |ATS03| . T> is row sparse). 

We want to find a quantum circuit which simulates J7[Y](i) = e~ lA ^ YS>t (note 
that because Y is regular the action of the adjacency matrix and laplacian are 
equivalent). To do this we apply the Trotter formula (see NCOO for details and 
further discussion) 

lim ( e - iAt / n e ~ iBt / n ^ n — & -i{A+B)t^ ^44-) 
n — >og 

By taking 0{n) — poly(| \B\, log(|V(F)|)) we gain a good approximation to 
the time evolution of A + B for time 0(t) = poly(n) |NC00| . 
Applying the Trotter formula to l|43|l we find that 

U[Y](t) RJ U e - iA ( x N)®Im.t/n\ e -iT> tiny, ^ Qr \J\Y\(t) ss ^ e ~ iA i x N)® J m t/n\ e -i"Dt/n\n _ 

(45) 

Because T> is row sparse, the simulation algorithm of ATS03 can be applied 
to simulate e~ lVt ^ n efficiently. (In order to guarantee the applicability of the 
Trotter formula we assume that ||A(Y)| grows polylogarithmically with |Y(Y)|.) 

We recursively reapply this construction to A(Xj^) (assuming, at each step, 
that the pull-back of the edges Xj is trivial for all but a small number of edges.) 
until we have expressed all instances of A(Y) with A(Xx). This construction fur- 
nishes an efficient (i.e. using 0(poly(log(|l / (Y)|))) elementary quantum gates) 
quantum circuit which simulates the propagator ?7[Y](t) accurately for times t 
which are 0(poly(log(|Y(Y)|))). 

The calculations in this section should be seen as representative of a general 
theory analogous to that initiated for quantum fourier transforms in [MR.R03 . 
Inat is, given a tower ot covering spaces Y > Ajy ► Ajv_i > ■ ■ ■ — > 



X%, what are the conditions the Xj must satisfy in order to give rise to an 
efficient gate sequence? We believe that such a theory is worth developing 
because it will potentially provide a class of quantum circuits whose behaviour 
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may be interesting from an algorithmic point of view. Interestingly, except 
for Cayley graphs, such quantum circuits will be unrelated to discrete fourier 
transforms. 



7 The Hidden Cover Problem 

The discussion in the previous section indicates that continuous-time quantum 
walks on certain graphs admit efficient gate decompositions whose size doesn't 
depend on the length of time the walk is simulated for. This feature can be 
exploited to give polynomial time (in log(|V(y)|)) quantum algorithms which 
can measure the hamiltonian A(Y). As a simple corollary of this we obtain an 
alternative observable for the hidden subgroup problem which is not immedi- 
ately equivalent to that employed by Shor's algorithm (and related quantum 
algorithms). Motivated by these new (efficiently implementable) observables 
we propose another generalisation of the hidden subgroup problem: the hidden 
covering space problem HiddenCover. 

The (coherent sampling) hidden subgroup problem for a group G (see, for 
example, NCOOj and references therein, for a discussion of the hidden subgroup 
problem) consists of a black box which outputs at random quantum states \i/) g H) 
which are equal superpositions of elements of cosets of a hidden subgroup H < 
G. (The Hilbert space H in this case is taken to be the group algebra C[G].) 
The problem is to determine H using as few of the coset states \ipgii) and as 
little quantum computational time as possible. 

The solution of the specific case G = "L/pqTi, where p and q are prime and 
H = Z/pZ or Z/gZ, is well-known — this is Shor's factoring algorithm. 

In the following we refer to a quantum state \<f>) in C[Z/nZ], n — pq, as a 
constant- coset state on cosets of a subgroup Z/gZ < Z/pqZ if it can be written 

p-i 

M=J>l«i>> (46) 

3=0 

where \otj) = l/y^ELo \j + l P)> 3 = 1>- •• >P~ 1 - 

The expansion of a constant coset state \<p) in the basis | W(j)) can be found 
via (this is essentially the discrete Poisson summation formula Tcr99 ): 



e p<7 j ■ — 



J2e-^ lk (47) 



fa,\qi A = 0, 



~ — 3k 



The expansion coefficients (W(k)\aj) are nonzero only when k is a multiple of 
9- 
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In the abelian hidden subgroup problem we have a black box which outputs 
at random the constant-coset states \ctj). The standard solution of the hidden 
subgroup problem proceeds by applying the quantum fourier transform to the 
\a,j) and then measuring in the computational basis. This yields an approx- 
imation r/q to the number r/q, for random integer r. After enough samples 
the identity of the hidden subgroup can be inferred by applying the continued- 
fractions algorithm to r/q. 

We now supply an alternative procedure to the standard quantum fourier 
transform which uses a quantum walk on the Cayley graph X (Z/pgZ, {±1}). 
We don't claim that this is any different to the standard quantum fourier trans- 
form algorithm for the HSP on Z/pgZ. The point is that the generalisation of 
this procedure to other graphs will not be equivalent to the quantum fourier 
transform method. At the moment, however, we can only perform the algorithm 
for cyclic groups. 

The eigenvalues of the hamiltonian (i.e. the Laplacian A) for the quantum 
walk on the cycle X(Z/pq'Z, {±1}) are given by 

X i = cos (^i) • («) 

Consequently, if the hamiltonian A [X] is measured exactly on a constant-coset 
state \4>) then, by the discussion surrounding (|4*?|) . the only eigenvalues that 

can be measured are those of the form Xj q = cos^^j^, for some random 

0<3 <q~l. 

We can effectively measure A[X] using J7[X](i) — which, as discussed in 
SjHl can be implemented efficiently, using the discretisation of von N eumann's 
prescription for measuring a hermitian operator given by |CDF+02| . Imple- 
menting this measurement yields an approximation Xj to an eigenvalue of A[X] 
for random j. Because cos(x) is continuous, the function 

i7^ = cos- 1 (A J ) (49) 

is a good approximation to the ratio j /q, for random < j < q — 1. Applying 
the continued fractions algorithm yields q. 

The previous result is suggestive of generalisations in the following way. 

Imagine we have a graph Y which is a covering space tt : Y — ► X for X. 
Imagine, further, we have a black box which outputs at random fibre-constant 
states \iP-k)- Recall that any fibre-constant state can be written in terms of the 
pull-backs of the eigenstates of X: 

\V(X)\-1 

\^)= E ^\E 3 {X)), (50) 

3=0 

where \Ej(X)) are the eigenstates of X. 

This means that if the hamiltonian A[Y] is measured on then it can only 
report eigenvalues of X, rather than eigenvalues from the full possible spectrum 
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of Y . If the spectrum of X(X) can be distinguished from the spectra X(Zi) of 
the other graphs Zi that Y is a covering space for then this measurement can 
identify the hidden graph X . As we showed previously, in the case where Y = 
X(Z/pqZ,{±l}), and X = X(Z/qZ,{±l}), then this is equivalent to solving 
the abelian hidden subgroup problem. 

Thus, generalising boldly, we believe that quantum computers ought to be 
able to efficiently solve the following problem. 



The hidden covering space problem Hidden Cover 
Input: 

f . A class C of mathematical objects. 

2. A quantisation scheme which associates a unitary matrix U(Y), with 
each Y e C, acting on an associated Hilbert space H(Y). 

3. An object Y from C with the promise that Y is a covering space 
tti : Y — > Xi only for a (known) set of objects Xi, I = 0, . . . , n. 

4. A black box which randomly emits fibre-constant states for some 
(unknown) projection 717, for some I = 0, . . . , n. 

Task: Determine the projection 717, and hence identify the base space Xi. 



We are willing to conjecture that HiddenCover is efficiently solvable on a 
quantum computer for certain classes of simple graphs using continuous-time 
quantum walks. 

A natural question arises: what happens when we apply the algorithms 
sketched above to the Cay ley graph for the dihedral group D n — (s,t\s n = 
t 2 = e)? In this case, we are after a hidden transposition (e, sH). Unfortunately, 
it can be verified that the Schreier graphs for the n hidden transpositions are 
isospectral, so that, with one fibre-constant state \ip n ) it is impossible to tell 
which transposition (e, sH) is hidden. 

8 Conclusions and Future Directions 

In this paper we have explored two related ideas, both of which are connected 
with covering space structures. The first is that a quantum evolution for an 
object Y which is a covering space for another object X can be equivalent to an 
evolution on the smaller object X. The second is that covering space structures 
can be exploited to give efficient gate sequences for their quantisations. The 
first idea can be exploited to give hitting-time speedups for such evolutions as 
quantum walks. Additionally, the first idea leads to the notion that hidden 
covering spaces can be identified spectrally. The second leads to the notion that 
quantum computers can do this efficiently 

We shall conclude with a list of future directions. 
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1. For the dihedral group consider walking on a graph Y which is not a 
simple cartesian product of log|D n | copies of the Cay ley graph of D n 
(this is what happens when you measure the propagator log|D n | times). 
Possible graphs to try might be certain graph products of the Cayley graph 
X(D n , S) for the dihedral group with generating set S. 

2. What about other mathematical objects? There are some promising can- 
didates, such as algebraic number fields, smooth manifolds, and knots 
which have natural structures amenable to quantisation. 

3. What sorts of computational problems are expressible as variants of Hid- 
denCover? Are there any interesting computational problems? 

Finally, a word on discrete-time quantum walks. Discrete-time quantum 
walks represent another quantisation scheme for simple graphs where the topol- 
ogy of the graph is encoded in the unitary quantisation. It is natural to ask 
how the ideas of this paper extend to the discrete-time quantum walks? It is 
interesting to remark that a coined quantum walk on a graph X is in fact a 
discrete-time quantum walk on a certain directed graph Y (namely Y is the 
line digraph of X — see, for example, GKWS98] for an application of line 
digraph in the design of algorithms) which is homomorphic to X . It can be 
observed that Y covers X (the spectrum of Y is the spectrum of X plus a zero 
eigenvalues with the appropriate multiplicity |Ros01p . Now, the construction 
mentioned above is only one possible quantisation of graphs. Are there other 
natural quantisations? Perhaps one can pick out a canonical quantisation by 
demanding that it respect covering space structures? 

Given a matrix M (over any field), a directed graph X is said to be the 
graph of M if the uv-th entry of M is nonzero if and only if there is a directed 
arc (u, v) for every pair of vertices u, v. A characterisation of graphs of uni- 
tary mat rices is still missing (see, for example, the open problems section of 
| ZKSS03p . In general, what is the relation between graphs of unitary matrices 
and covering spaces? The study of this problem may be useful in understand- 
ing the combinatorial properties of certain combinatorial designs like weighing 
matrices. 
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